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While the baryon asymmetry of the Universe is nowadays well measured by cosmological obser- 
vations, the bounds on the lepton asymmetry in the form of neutrinos are still significantly weaker. 
We place limits on the relic neutrino asymmetries using some of the latest cosmological data, taking 
into account the effect of flavor oscillations. We present our results for two diff'erent values of the 
neutrino mixing angle and show that for large 613 the limits on the total neutrino asymmetry 
become more stringent, diluting even large initial flavor asymmetries. In particular, we find that the 
present bounds are still dominated by the limits coming from Big Bang Nucleosynthesis, while the 
limits on the total neutrino mass from cosmological data are essentially independent of 613. Finally, 
we perform a forecast for COrE, taken as an example of a future CMB experiment, and find that it 
could improve the limits on the total lepton asymmetry approximately by up to a factor 6.6. 

PACS numbers: 98.80.-k, 14.60. Pq. 26.35. +c, 98.70.Vc, 98.80.Es 



I. INTRODUCTION 

Quantifying the asymmetry between matter and anti- 
matter of the Universe is crucial for understanding some 
of the particle physics processes that might have taken 
place in the early Universe, at energies much larger than 
the ones that can be reached currently in particle acceler- 
ators. Probes of the anisotropics of the cosmic microwave 
background (CMB) together with other cosmological ob- 
servations have measured the cosmological baryon asym- 
metry rjb to the percent level thanks to very precise mea- 
surements of the baryon density [l| . For the lepton asym- 
metries, while they are expected to be of the same order 
of the baryonic one due to sphaleron effects that equili- 
brate both asymmetries, it could be the case that other 
physical processes lead instead to leptonic asymmetries 
much larger than r^, (see, e.g., [IHSDj with consequences 
for the early Universe phase transitions Q , cosmological 
magnetic fields Q , and the dark matter relic density 0- 
[^. Neutrino asymmetries are also bound to be nonzero 
in the presence of neutrino isocurvature perturbations, 
like those generated by curvaton decay |10l - [l3 |. Those 
large neutrino asymmetries could have been imprinted in 
the cosmological data [ij, [ij] , and although the limits 
on such asymmetries have been improving over the last 
years, current constraints are still many orders of magni- 
tude weaker than the baryonic measurement. 

On the other hand, thanks to the neutrino oscilla- 
tions the initial primordial flavor asymmetries are redis- 
tributed among the active neutrinos before the onset of 



Big Bang Nucleosynthesis (BBN) [Tl-[i3> 

which makes 

the knowledge of the oscillation parameters important 
for correctly interpreting the limits on such asymme- 
tries. Nowadays all of those parameters are accurately 
measured (see e.g. [H, [l^), with the exception of the 
mixing angle ^13 that only recently started to be signifi- 
cantly constrained. In fact, several neutrino experiments 
over the last year gave indications of nonzero values for 
sin^ 6*13 |20l - [23 | , and recently the Daya Bay reactor ex- 
periment claimed a measurement of sin^(20i3) = 0.092 ± 
0.016(stat.)±0.005(syst.) at 68% C.L. [2I, excluding a 
zero value for ^13 with high significance. The same find- 
ing has been also reported by the RENO Collaboration 
H, sin2(26li3) = 0.113 ±0.013(stat.)±0.019(syst.) (68% 



C.L.). 

Finally, yet another important piece of information for 
reconstructing the neutrino asymmetries in the Universe 
is the measured value of the relativistic degrees of free- 
dom in the early universe, quantified in the so-called 
effective number of neutrinos, iVoff. In the case of the 
three active neutrino flavors with zero asymmetries and 
a standard thermal history, its value is the well-known 
Ncfi — 3.046 [2^, but the presence of neutrino asymme- 
tries can increase that number while still satisfying the 
BBN constraints [2^ . Interestingly enough, recent CMB 
data has consistently given indications of A^off higher 
than the standard value: recently the Atacama Cosmol- 
ogy Telescope (ACT) [Hi and the South Pole Telescope 
(SPT) [H, [li] have found evidence for N^s > 3.046 at 
95% C. L., making the case for extra relativistic degrees 
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of freedom stronger (see also [30| ) • It should however be 
kept in mind that other physical processes, like e.g. the 
contribution from the energy density of sterile neutrinos 
[31I [32! or of gravitational waves [SJ] , could also lead to 
a larger value for A^off ■ 

Some recent papers have analyzed the impact of neu- 
trino asymmetries with oscillations on BBN [2^ [13, , 
mainly because data on light element abundances domi- 
nate the current limits on the asymmetries. Some studies 
using CMB data can be found in the literature (see for 
instance [36l - l38j for limits on the degeneracy parameters 
^1, using the WMAP data and ^ for the effect of the 
primordial Helium fraction in a Planck forecast), but our 
paper improves on that in two directions. First, we used 
for our analysis the neutrino spectra in the presence of 
asymmetries after taking into account the effect of flavor 
oscillations. Second, we checked the robustness of our 
results comparing the analysis of CMB and BBN data 
with a more complete set of cosmological data, including 
in particular supernovae la (SNIa) data (i^l, the mea- 
surement of the Hubble constant from the Hubble Space 
Telescope (HST) [U, and the Sloan Digital Skv Survey 
(SDSS) data on the matter power spectrum |42|. While 
current CMB measurements and the other datascts are 
not expected to improve significantly the constraints on 
the asymmetries, they constrain the sum of the neutrino 
masses, giving a more robust and general picture of the 
cosmological parameters. 

Our goals in this work is twofold: first, we constrain 
the neutrino asymmetries and the sum of neutrino masses 
for both zero and nonzero values of ^13 using some of the 
latest cosmological data to obtain an updated and clear 
idea of the limits on them using current data; second, 
we perform a forecast of the constraints that could be 
achievable with future CMB experime nts, taking as an 
example the proposed mission COrE^ [4^. Given that 
current constraints are basically dominated by the BBN 
constraints, we use our forecast to answer the more gen- 
eral question of whether future CMB experiments can be 
competitive with the BBN bounds. 

This paper is organized as follows. Initially, we briefly 
review in Sec. |TT] the dynamics of the neutrino asymme- 
tries prior to the BBN epoch. With those tools in hand, 
we proceed to study in Sec. IHII the impact on cosmo- 
logical observablcs of the neutrino asymmetries for two 
values of the mixing angle ^13 using current cosmologi- 
cal data. We then step towards the future and describe 
in Sec. IIVI our forecast for the experiment COrE, where 
we study the potential of the future data from Icnsing of 
CMB anisotropics to constrain some of the cosmological 
parameters (in particular, neutrino asymmetries and the 
sum of the neutrino masses) with great precision. Finally, 
in Sec. |V]we draw our conclusions. 



II. EVOLUTION OF COSMOLOGICAL 
NEUTRINOS WITH FLAVOR ASYMMETRIES 

The dynamics of the neutrino distribution functions in 
the presence of flavor asymmetries and neutrino oscilla- 
tions in the earl y U niverse has been discussed in detail 
in the literature p6l . Is^ . HH , and here we will only briefly 
review its main features and its consequences for the late 
cosmology. 

We assume that flavor neutrino asymmetries, rj^^ , were 
produced in the early Universe. At large temperatures 
frequent weak interactions keep neutrinos in equilibrium 
thus, their energy spectrum follows a Fermi-Dirac distri- 
bution with a chemical potential fXi,^ for each neutrino 
fiavor. If = P-Vc/T is the degeneracy parameter, the 
asymmetry is given by 
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12C(3) 



(1) 



Here n^^ ('t-Po) denotes the neutrino (antineutrino) num- 
ber density, is the photon number density, and C(3) = 
1.20206. 

As usual, we will write the radiation energy density of 
the Universe in terms of the parameter A^eff 1 the effective 
number of neutrinos, as 



Pr = P-1 



A/?. 



N,- 



(2) 



with Ncs ~ 3.046 the value in the standard case with zero 
asymmetries and no extra relativistic degrees of freedom 
[25|. Assuming that equilibrium holds for the neutrino 
distribution functions, the presence of fiavor asymmetries 
leads to an enhancement 
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^ http://www.corc-mission.org 



Note that a neutrino degeneracy parameter of order 
^ 0.3 is needed in order to have a value of AA'oft at 
least at the same level of the effect of non-thermal distor- 
tions discussed in [2^. This corresponds to 77^^ ~ 0(0.1). 
On the other hand, the primordial abundance of ^He de- 
pends on the presence of an electron neutrino asymmetry 
and sets a stringent BBN bound on rj^^ which docs not 
apply to the other flavors, leaving a total neutrino asym- 
metry of order unity unconstrained jH, |45|. However, 
this conclusion relies on the absence of effective neu- 
trino oscillations that would modify the distribution of 
the asymmetries among the different flavors before BBN. 

The evolution of the neutrino asymmetries in the 
epoch before BBN with three-flavor neutrino oscillations 
is found by solving the equations of motion for 3 x 3 
density matrices of the flavor neutrinos as described in 
[46l [43, including time-dependent vacuum and matter 
terms, both from background and neutrinos, as well 
as the collision integrals from neutrino weak interactions. 
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TABLE I: Cosmological and neutrino parameters. 



Type 


Symbol 


Meaning 


Uniform Prior 


Primary 




Baryon density 


(0.005,0.1) 


Cosmological 




Dark matter density" 


(0.01,0.99) 


Parameters 


r 


Optical depth to reionization 


(0.01,0.8) 






Angular scale of the sound horizon at the last scattering 


(0.5, 10) 




ris 


Scalar index of the power spectrum 


(0.5,1.5) 




log [10"' A«] 


Scalar amplitude of the power spectrum ' 


(2.7,4) 


Neutrino 


mi(eV) 


Mass of the lightest neutrino 


(0,1) 


Parameters 




Total asymmetry at T = 10 MeV 


(-0.8,0.8) 






Initial electron neutrino asymmetry at T = 10 MeV 


(-1.2,1.2) 


Derived 


h 


Reduced Hubble constant'' 




Parameters 


ATVeff 


Enhancement to the standard effective number of neutrinos'^ 





"Also includes neutrinos. 

''at the pivot wavenumbcr kg = 0.05 Mpc~^. 
'^We assume here normal hierarchy. 
rf_f/o = lOO/i km s-i Mpc-i. 
^AT^fj = 3.046. 



This was done under certain approximations in refs. [isl - 
[iTj . where it was shown that neutrino oscillations are 
indeed effective before the onset of BBN. Therefore, the 
total lepton asymmetry is redistributed among the neu- 
trino flavors and the BBN bound on r^^^ can be translated 
into a limit on iji, = r/^^ + r/^^ + t]^^ , unchanged by oscilla- 
tions and constant until electron-positron annihilations, 
when it decreases due to the increase in the photon num- 
ber density. 

The temperature at which flavor oscillations become 
effective is important not only to establish r/^^ at the 
onset of BBN, but also to determine whether weak in- 
teractions with e^e~ can still keep neutrinos in good 
thermal contact with the primeval plasma. Oscillations 
redistribute the asymmetries among the flavors, but only 
if they occur early enough interactions would preserve 
Fermi-Dirac spectra for neutrinos, in such a way that the 
degeneracies are well defined for each rj^^ and the re- 
lation in Eq. ([3]) remains valid. This the case of early 
conversions of muon and tau neutrinos, since oscillations 
and collisions rapidly equilibrate their asymmetries at 
T ~ 15 MeV [l5|. Therefore one can assume the initial 
values 7y™ = ry™ = 77™, leaving as free parameters 77"' 
and the total asymmetry rji, = 77™ -I- 2?7™ . 

If the initial values of the flavor asymmetries jyjf and 
?7™ have opposite signs, neutrino conversions will tend 
to reduce the asymmetries which in turn will decrease 
iVeff. But if flavor oscillations take place at temperatures 
close to neutrino decoupling this would not hold and an 
extra contribution of neutrinos to radiation is expected 
with respect to the value in Eq. ([3]), as emphasized in 
(26j and shown in Fig. [1] where the A^eff isocontours for 
non-zero mixing are compared with those obtained from 
the frozen neutrino distributions taking into account the 
effect of flavor oscillations [s^. One can see that oscil- 
lations efficiently reduce A'cs for neutrino asymmetries 



with respect to the initial values from Eq. 

The evolution of the neutrino and antincutrino distri- 
bution functions with non-zero initial asymmetries, from 
T = 10 MeV untn BBN, has been calculated in Hi,!!!]. 
Here we use the final numerical results for these spectra 
in a range of values for 77™ and rji, as an input for our anal- 
ysis, described in the next Section. Note that an analysis 
in terms of the degeneracy parameters ^q, as done for in- 
stance in [s^ is no longer possible. We adopt the best fit 
values for the neutrino oscillation parameters quoted in 
[m , assuming a normal hierarchy of the neutrino masses, 
except for the mixing angle ^13, for which we will adopt 
two distinct values: 6'i3 = and sin^ 6'i3 = 0.04. The lat- 
ter is close to the upper limit placed by the Daya Bay [2^ 
and RENO [13] experiments on this mixing angle (with 
a best-fit value of sin^ 6I13 = 0.024 and sin^ 6*13 = 0.029, 
respectively), and is used as an example to understand 
the cosmological implications of a nonzero 6*13. Moreover, 
since the flavor asymmetries equilibrate for large values 
of this mixing angle, the cosmological effects are similar 
for sin^ 6*13 > 0.02, as in the case of an inverted hierarchy 
for a broad range of 6*13 values (see, for instance. Fig. 4 
of Rcf. [13). As for the case ^13 = 0, though it seems 
presently disfavoured with a high statistical significance 
after the Daya Bay and RENO results, we have decided 
to include it for comparison. 



III. COSMOLOGICAL CONSTRAINTS ON 
NEUTRINO PARAMETERS 

Having set the basic framework for the calculation of 
the neutrino distribution functions in the presence of 
asymmetries and for different ^13, we can now proceed 
to investigate its cosmological effects. 

In order to constrain the values of the cosmological 
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FIG. 1: Final contribution of neutrinos with primordial asymmetries to the radiation energy density. The isocontours of NcB 
on the plane 77"^ vs. t)^, including flavor oscillations, are shown for two values of sin^ ^13: (blue solid curves, left panel) and 
0.04 (red solid curves, right panel) and compared to the case with zero mixing (dashed curves). The dotted line corresponds 
to rjv — ri,j^ {x — h,t), where one expects oscillations to have negligible effects. 



neutrino asymmetries, we compare our results to the ob- 
servational data. In particular, we use a modified ver- 
sion of the CAMB code^ [13 to evolve the cosmological 
perturbations and obtain the CMB and matter power 
spectra in the presence of non-zero neutrino asymmetries 
in the neutrino distribution functions. We checked that 
the spectra computed by our modified CAMB version 
are consistent up to high accuracy with those obtained 
with CLASS [43, that incorporates the models consid- 
ered here in its public version. This version of CAMB 
is interfaced with the Markov chain Monte Carlo pack- 
age CosmoMC^ [lOl that we use to sample the parameter 
space and obtain the posterior distributions for the pa- 
rameters of interest. 

Wc derive our constraints in the framework of a fiat 
ACDM model with the three standard model neutrinos 
and purely adiabatic initial conditions. The parameters 
we use are described in Table U as well as the range of 
the flat priors used. As can be seen, six of them are 
the standard ACDM cosmological parameters, and we 
add to those three new parameters, namely the mass of 
the lightest neutrino mass cigcnstatc mi (the other two 
masses are calculated using the best fit for Am|i and 
Am§i obtained in [isj . assuming normal hierarchy) and 
the two neutrino asymmetries we mentioned earlier, yyjf 
and 771/. The values of the effective degeneracy parame- 
ters £,a after BBN'', needed by CAMB, are pre-calculated 



as a function of the asymmetries (following the method 
described in the previous section) over a grid in (77™ , rj^) 
and stored on a table, used for interpolation during the 
Monte Carlo run. 

A comment on the parameterization is in order. It 
is a standard practice in cosmological analyses to pa- 
rameterize the neutrino masses via Q,^h? or equivalently 
fi, = D.^/D.dm, and from that (assuming that neutri- 
nos decoupled at equilibrium) derive the sum of neutrino 
masses, which are taken to be degenerate. The presence 
of lepton asymmetries dramatically changes this simple 
scheme. Now the neutrino number density is a compli- 
cated function of the ?7's obtained from a non-equilibrium 
distribution function. When fi, is used, any effect related 
to the way in which the total neutrino density is shared 
among the different mass eigenstates is completely lost. 
In that sense, the parameterization used in this paper 
looks more physically motivated since energy densities of 
neutrinos arc constructed from two fundamental quan- 
tities, namely their phase space distributions and their 
masses. 

The most basic dataset that we consider only con- 
sists of the WMAP 7-year temperature and polarization 
anisotropy data. We will refer to it simply as "WMAP" . 
The likelihood is computed using the the WMAP like- 
lihood code publicly available at the LAMBDA web- 



http : //camb. Inf o/| Fermi-Dirac-like functions with an effective and temperature 



http : //cosmologist . inf o/cosmomc/ 1 Ta [3J| , which are related to the first two moments of the distri- 



The neutrino distribution functions can be parameterized by bution, the number density and energy density. 
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site^. We marginalize over the amplitude of the Sunyaev- 
Zel'dovich signal. 

In addition to the WMAP data, we also include the 
BBN measurement of the ^He mass fraction Yp from the 
data collection analysis done in jsij . in the form of a 
Gaussian prior 

yp = 0.250 ±0.003 (Icr). (4) 

Indeed, some authors have recently reported a larger 
central value, Yp ^ 0.257 f52| - [5^. with quite differ- 
ent uncertainty determinations. In [55| using a Markov 
chain Monte Carlo technique already exploited in [s^ ]. 
the primordial value of ''He decreased again to Yp = 
0.2534 ±0.0083, which is compatible at Icr with Q. We 
will not use these results in our analysis, but we will com- 
ment on their possible impact in the following. We also 
note that in [5a] a robust upper bound Yp < 0.2631 (95 % 
C.L.) has been derived based on very weak assumptions 
on the astrophysical determination of ^He abundance, 
namely that the minimum effect of star processing is 
to keep constant the helium content of a low-metallicity 
gas, rather than increase it, as expected. As we will 
show, the measurement of Yp currently dominates the 
constraints on the asymmetries: if we were to conserva- 
tively allow for larger uncertainties on that measurement, 
like for example those reported in (55| . our constraints 
from present data would correspondingly be weakened. 
Moreover, we decided not to use the Deuterium mea- 
surements since at the moment they are not competitive 
with Helium for constraining the asymmetries (see, e.g., 
Fig. 6 of Ref. [IJl), although there are recent claims that 
they could place strong constraints on Nes at the level of 
AA'off — ±0.5 [53|. This is a very interesting perspective 
but at the moment, Deuterium measurements in different 
QSO absorption line systems show a significant disper- 
sion, much larger than the quoted errors. 

The dataset that uses both WMAP 7-year data and the 
determination of the primordial abundance of Helium as 
in dH) will be referred to as "WMAP+He". Measure- 
ments of Yp represent the best "leptometer" currently 
available, in the sense that they place the most stringent 
constraints on lepton asymmetries for a given baryonic 
density [11]. The ''He mass fraction depends on the bary- 
onic density, the electron neutrino degeneracy parameter 
and the effective number of neutrino families. Thus, in 
order to consistently implement the above determination 
of Yp in our Monte Carlo analysis, we compute AA'off 
and coming from the distribution functions calculated 
with the asymmetries (as explained in the previous sec- 
tion) and store them on a table. During the CosmoMC 
run, we use this table to obtain by interpolation the val- 
ues AA'off and corresponding to given values of the 
asymmetries (which are the parameters actually used in 
the Monte Carlo) , and finally to obtain Yp as a function of 
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^Ncff, and Vt^h?'. Notice that this approach is slightly 
less precise than the one used in Refs. [H, [s^, where 
a full BBN analysis was performed, but this approxima- 
tion should suffice for our purposes, especially taking into 
account that we will be comparing BBN limits on the 
asymmetries with the ones placed by other cosmological 
data, that as we shall see are far less constraining. In 
any case, we have checked that the agreement between 
the interpolation scheme and the full BBN analysis is at 
the percent level. 

We derive our constraints from parallel chains gener- 
ated using the Metropolis-Hastings algorithm. For a sub- 
set of the models, we have also generated chains using the 
slice sampling method, in order to test the robustness of 
our results against a change in the algorithm. We use 
the Gelman and Rubin R parameter to evaluate the con- 
vergence of the chains, demanding that i? — 1 < 0.03. 
The one- and two-dimensional posteriors are derived by 
marginalizing over the other parameters. 

Our results for the cosmological and neutrino param- 
eters from the analysis are shown in Table [Hi while 
Fig. [2] shows the marginalized one-dimensional probabil- 
ity distributions for the lightest neutrino mass, the initial 
electron-neutrino asymmetry, and the total asymmetry, 
for the different values of 6'i3. Notice that the posterior 
for (middle panel) is still quite large at the edges of 
the prior range. This happens also for both the 77™ and 
■q^ posteriors obtained using only the WMAP data (not 
shown in the figure). Since the priors on these parame- 
ters do not represent a real physical constraint (as in the 
case TOj. > 0), but just a choice of the range to explore, 
we refrain from quoting 95% credible intervals in these 
cases, as in order to do this one would need knowledge of 
the posterior in all the region where it significantly differs 
from zero. However, it is certain that the actual 95% C.I. 
includes the one that one would obtain using just part of 
the posterior (as long as this contains the peak of the dis- 
tribution). If we do this, we obtain constraints that are 
anyway much worse than those from BBN. Finally, we 
also stress that if a larger experimental determination of 
Yp or measurements with larger uncertainities were used, 
as those reported in (s^ - ls^ . BBN would show a prefer- 
ence for larger values of A^eff as well. 

Concerning the neutrino asymmetries, shown in the 
middle and right panels of Fig. [21 we notice that while the 
initial flavor asymmetries remain highly unconstrained 
by current data, the total asymmetry constraint improves 
significantly for 6*13 ^ 0. This result agrees with previ- 
ous results from BBN-only studies [13, [1^ , and it is a 
result of the equilibration of flavor asymmetries when 
6'i3 is large (see, e.g., Fig. 5 of Ref. [13]). When the fla- 
vors equilibrate in the presence of a nonzero mixing angle 
(sin^ 013 = 0.04 in our example) the total asymmetry is 
distributed almost equally among the different flavors, 
leading to a final asymmetry 77^" « 77^" « -q^/?) (where 
X = ijl,t). Hence, the fact that the BBN prior requires 
77^" ~ for the correct abundance of primordial Helium 
(see Fig. [3]) leads to a strong constraint on the constant 
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FIG. 2: One-dimensional posterior probability density for mi, rj]^^, and 77,^ for the WMAP+He dataset. 



TABLE II: 95% C.L. constraints on cosmological parameters for the WMAP and WMAP+He datasets. 



Parameter 


WMAP 




WMAP±He 






sin^ 6*13 = 


sin^ 6*13 = 0.04 


sin^ 6I13 = 


sin^ 6*13 = 0.04 


100 flbh'^ 


2 20+°'^* 


2 20+° " 


2.20 ±0.12 


2.20 ±0.12 


r 

woe. 

Us 

log \W^°A,] 


0.118 ±0.016 

0.085tH^^ 
1.0387 ± 0.0063 

0.953 ± 0.032 


0.085t°:«f, 

±.uoot)_o 0063 
U.yOO_o.o33 

3 062+° °*° 

O.UU^_0 079 


0.119 ±0.017 
0.085l°:°i? 

l-038ll°:°S 
0.955i°:°lt 
3.0681°:°?^ 


0. 117 ±0.016 
0.085^0 027 

1. U.3»C_0.0054 

Qt^2+° °3i 

U.»O^_0.032 

3 062+° °''^ 

O.UUZ_0 075 


mi (eV) 


< 0.39 

a 


< 0.38 

a 


< 0.38 

a 


< 0.38 

a 


< 


a 


a 


[-0.64; 0.72] 


[-0.071; 0.054] 


h 

AN,B 


U.DOZ_g Qg3 

< 0.32 


fiS3+° °*^ 
< 0.16 


656+° °*^"' 

U.DOD„o.o81 

< 0.43 


0.6501°:°^? 
< 0.03 



"The 95% confidence region is not well-defined in these cases be- 
cause the posterior does not vanish at the end of the prior range 
(see e.g. the middle panel of Fig. [2]l. See discussion in the text. 



total asymmetry, -0.071 < rj^ < 0.054 (95% C.L.). 

On the other hand, since the constraints come most 
from the distortion in the electron neutrino distribution 
function, when ^13 = (and therefore there is less mix- 
ing) the direct relation between 77^" and rj^ is lost. In 
this case, the total asymmetry could still be large, even 
if the final electron neutrino asymmetry is small, as sig- 
nificantly asymmetries can still be stored on the other 
two flavors, leading to a constraint an order of magni- 
tude weaker than the previous case, —0.64 < rjy < 0.72 
(95% C.L.). As expected, this is reflected on the allowed 
ranges for AA'off , as shown in Fig.Hl while for 6^13 = the 
AiVeff ~ 0.5 are still allowed by the data, nonzero values 
of this mixing angle reduce the allowed region in the pa- 
rameter space by approximately an order of magnitude 
in both AA^off and 77^. 

We confirmed in our analysis that the constraints on 
the asymmetry are largely dominated by the BBN prior 
at present. This is shown in Fig. [5l where we compare 
the results of our analysis with a more complete dataset 
(which we refer to as ALL) that includes distance mea- 
surements of SNIa from the SDSS compilation |40| and 
the HST determination of the Hubble constant Hq [4l| , as 
well as data on the power spectrum of the matter density 



field, as reconstructed from a sample of Luminous Red 
Galaxies of the SDSS Seventh Data Release [42|. This 



is due to the fact that other cosmological data constrain 
the asymmetries via their effect on increasing ./Voff, and 
currently the errors on the measurement of the effective 
number of neutrinos [l|, [27l - l29j are significantly weaker 
than our prior on Yp, eq. (|4])^. The fact that bounds on 
Icptonic asymmetries are dominated by the BBN prior 
(i.e. by ''He data) is also confirmed by the similarity of 
our bounds on (77^, 77™) with those of [1^. Note that 
the limits reported in [35] sound weaker, because they 
are frequentist bounds obtained by cutting the parame- 



" On the other hand, these other cosmological data sets have an 
impact on other parameters like e.g. the neutrino mass. But 
since in this work we are primarily interested in bounding the 
asymmetries, we prefer to stick to the robust WMAP-|-He data 
set. In that way, our results are not contaminated by possible 
systematic uncertainties in the other data. Actually, the inclu- 
sion of all external datasets (in particular, of SNIa together with 
Hq) reveals a conflict between them, leading to a bimodal pos- 
terior probability for Q.ii^h? and to a preference for mi > at 
95% C.L. 
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FIG. 3: 68% and 95% confidence regions in total neutrino 
asymmetry rj^, vs. the primordial abundance of Helium Yp 
plane for On = (blue) and sin^ On = 0.04 (red), from the 
analysis of the WMAP-|-He dataset. Notice the much stronger 
constraint for the nonzero mixing angle due to the faster equi- 
libration of flavor asymmetries. 




FIG. 4: Two-dimensional 68% and 95% confidence regions 
in the {r]v, Ncs) plane from the analysis of the WMAP-|-He 
dataset, for ^13 = (blue) and sin^ 6^13 — 0.04 (red). Even 
for zero S13 the data seem to favor A'oft around the standard 
value iVeff = 3.046. 



tcr probability at Ax^ = 6.18, i.e. they represent 95% 
bounds on joint two-dimensional parameter probabilities 
(in the Gaussian approximation). The one-dimensional 
95% confidence limits, corresponding to Ax^ = 4, are 
smaller and very close to the results of the present paper. 
We also checked that using our codes and data sets, we 
obtain very similar results when switching from Bayesian 
to frequentist confidence limits. 

We conclude this section noting that the current con- 
straints on the sum of neutrino masses are robust un- 
der a scenario with lepton asymmetries, as those extra 
degrees-of-freedom do not correlate with the neutrino 
mass. On the other hand, to go beyond the BBN limits on 
the asymmetries more precise measurements of A^eff are 
clearly needed, and in the next section we forecast the re- 
sults that could be achievable with such an improvement 
using COrE as an example of future CMB experiments. 



" 0.4 
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FIG. 5: One-dimensional posterior probability density for r\i, 
comparing the WMAP-t-He and the ALL datasets. As men- 
tioned in the text, the constraints on the total asymmetry do 
not improve significantly with the inclusion of other cosmolog- 
ical datasets, as they are mainly driven by the determination 
of the primordial Helium abundance. 



IV. FORECAST 

Given that the current constraints on the lepton asym- 
metries are dominated by their effect on the primordial 
production of light elements, one can ask whether future 
cosmological experiments can improve over the current 
limits imposed by BBN. With that goal in mind, we take 
as an example a proposed CMB experiment, COrE (Cos- 
mic Origins Explorer) [4^, designed to detect the pri- 
mordial gravitational waves and measure the CMB grav- 
itational lensing deflection power spectrum on all linear 
scales to the cosmic variance limit. The latter is of special 
interest for this work, as the CMB lensing is expected to 
probe with high sensitivity the absolute neutrino masses 
and iVoff [13. 
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We used the package FuturCMB'' in combination with 
CAMB and CosmoMC for producing mock CMB data, 
and fit it with a hkehhood based on the potential sensitiv- 
ity of COrE. We include, also in this case, the information 
coming from present measurements of the Helium frac- 
tion, encoded in the Gaussian prior ^ . We consider five 
of COrE's frequency channels, ranging from 105 to 225 
GHz, with the specifications given in [43l | and reported for 
convenience in Table IHIl and assume an observed frac- 
tion /sky = 0.65. We do not consider other channels as 
they are likely to be foreground dominated. We take a 
maximum multipole imax = 2500. In our analysis, we 
have assumed that the uncertainties associated to the 
beam and foregrounds have been properly modeled and 
removed, so that we can only consider the statistical un- 
certainties. Those are optimistic assumptions, as under 
realistic conditions systematic uncertainties will certainly 
play an important role. In that sense, our results rep- 
resent an illustration of what future CMB experiments 
could ideally achieve. 



TABLE IV: Fiducial values for the cosmological parameters 
for the COrE forecast. 



Frequency [GHz] 


fi'twhm [arcmin] 




ap [fiK] 


105 


10.0 


0.268 


0.463 


135 


7.8 


0.337 


0.583 


165 


6.4 


0.417 


0.720 


195 


5.4 


0.487 


0.841 


225 


4.7 


0.562 


0.972 



TABLE III: Experimental specifications for COrE For 
each channel, we list the channel frequency in GHz, the 
FWHM in arcminutes, the temperature (or) and polariza- 
tion ((7p) noise per pixel in /iK. 



We use CMB lensing information in the way described 
in [g^l, assuming that the CMB lensing potential spec- 
trum will be extracted from COrE maps with a quadratic 
estimator technique. 

For the forecast we adopt the fiducial values for the cos- 
mological parameters shown in Table IIVI for both cases 
of 6*13 discussed previously. The two sets of fiducial val- 
ues correspond to the best-fit models of the WMAP-|-He 
dataset for the two values of ^13. In the case of the neu- 
trino mass, since the likelihood is essentially flat between 
and 0.2 eV, we have chosen to take toi = 0.02 eV. This 
is below the expected sensitivity of COrE and should 
thus be essentialy equivalent to the case where the light- 
est neutrino is massless. 

The sensitivities on the neutrino parameters for COrE 
are shown in Fig.[n]for the two values of 6*13. As expected 
for the sum of the neutrino masses, the constraints are 
significantly better than the current ones, and could in 
principle start probing the minimal values guaranteed by 
fiavor oscillations QW,. Note that our forecast error for 



PcLrcLineter 


Fiducial Value 


Fiducial Value 




fsin^ i9th =0) 


(sin^ 013 = 0.04) 




0.0218 


0.0224 




0.121 


0.118 




0.0873 


0.0865 


h 


0.709 


0.705 


Us 


0.978 


0.968 


log [W^°As] 


3.12 


3.08 


mi (eV) 


0.02 


0.02 


















TABLE V: 95% confidence intervals for the neutrino param- 


eters with COrE. 






Parameter 


sin^ 013 = 


sin^ 6I13 = 0.04 


mi (eV) 


< 0.049 


< 0.048 




[-0.20; 0.20] 


[-0.25; 0.24] 




[-0.12; 0.09] 


[-0.048; 0.030] 
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nil differs slightly from the one presented in [43|, most 
probably because the forecasts in this reference are based 
on the Fisher matrix approximation. But our main goal 
in this section is to discuss how COrE observations will 
help improving the limits on the asymmetries discussed 
previously, that are basically dominated by the available 
measurements of the ^He abudance. The right panel of 
Fig. [5] shows the forecasted posterior probability distri- 
bution for rji, , and the marginalized constraints for it are 
listed in Table IVl for both values of 6113; in particular, the 
vertical lines of the right panel show the 95% C.L. limits 
obtained from the full BBN analysis of Ref. [IH . Com- 
paring the values from Tables [H] and |V] one can see that 
an experiment like COrE would improve current 95% lim- 
its on the total leptonic asymmetry by nearly a factor 6.6 
(6*13 = 0) and 1.6 (sin^ 6*13 = 0.04), competitive over the 
constraints from ''He abundance only. It should be noted 
that the error bars on the primordial abundances are very 
difficult to be reduced due to systematic errors on astro- 
physical measurements [5l| . and therefore it is feasible 
that CMB experiments will be an important tool in the 
future to improve the constraints on the asymmetries. 
Notice however that, since the CMB is insensitive to the 
sign of the 77's, BBN measurements will still be needed 
in order to break this degeneracy. 

Finally, in Fig. [7] we show the COrE sensitivity on the 
asymmetries in the plane r]^ vs. jyjf compared to the 
constraints of Sec. Illll obtained using current data and to 
the full BBN analysis of Ref. [1^ . Notice that in the case 
013 = the constraints of the previous section are quite 
less constraining than the ones coming from the full BBN 
analysis because we are not using deuterium data, known 
to be important to close the contours on the asymmetries 
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FIG. 6: One-dimensional probability distribution function for m\ and r]^ for COrE forecast. The middle panel shows that an 
experiment like COrE could start constrain the initial electron neutrino asymmetry. The vertical lines on the right panel show 
the current 95% C.L. limits obtained in the previous section. The errors on the asymmetries are improved by approximately a 
factor 6.6 or 1.6 for S13 = and sin^ ^13 = 0.04, respectively, compared to the results shown in Fig. m 
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FIG. 7: The 95% C.L. contours on the 77^ vs. ?7J," plane from our analysis with current data (WMAP+He dataset, black 
dotted) compared to the results of the BBN analysis of Ref. [3^1 (blue dashed) and with the COrE forecast (red solid). 



plane, especially for small values of 613 [3J]. Moreover, 
future CMB experiments have the potential to reduce the 
allowed region, dominating the errors in this analysis. 



In summary, an experiment like COrE is capable of 
improving the constraints on the Icpton asymmetries by 
up to a factor 6.6 on the total and/or flavor asymmetries 
depending on the value of the mixing angle ^13. In ad- 
dition to that, such an experiment would also constrain 
other cosmological parameters (in particular the sum of 
the neutrino masses) with significant precision, providing 
yet another step towards the goal of accurately measuring 
the properties of the Universe. 



V. CONCLUSIONS 

Understanding the physical processes that took place 
in the early Universe is a crucial ingredient for deci- 
phering the physics at energies that cannot be currently 
probed in terrestrial laboratories. In particular, since the 
origin of the matter-antimatter is still an open question 
in cosmology, it is important to keep an open mind for 
theories that predict large lepton asymmetries. In that 
case, constraining total and fiavor neutrino asymmetries 
using cosmological data is a way to test and constrain 
some of the possible particle physics scenarios at epochs 
earlier than the BBN. 

For that, we initially used current cosmological data 
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to constrain not only the asymmetries, but also to un- 
derstand the robustness of the cosmological parameters 
(and the limits on the sum of the neutrino masses) for 
two different values of the mixing angle ^13 to account 
for the evidences of a nonzero value for this angle. Our 
results confirm the fact that at present the limits on the 
cosmological Icpton asymmetries arc dominated by the 
abundance of primordial elements generated during the 
BBN, in particular the abundance of *He, currently the 
most sensitive "leptometer" available. 

However, future CMB experiments might be able to 
compete with BBN data in what concerns constraining 
lepton asymmetries, although BBN will always be needed 
in order to get information on the sign of the 77's. We 
took as an example the future CMB mission COrE, pro- 
posed to measure with unpreccdcnt precision the lensing 
of CMB anisotropics, and our results indicate that it has 
the potential to significantly improve over current con- 
straints while, at the same time placing limits on the 
sum of the neutrino masses that are of the order of the 
neutrino mass differences. 

Finally, we notice that for the values of 6*13 measured 
by the Daya Bay and RENO experiments the limits on 
the cosmological lepton asymmetries and on its asso- 
ciated effective number of neutrinos are quite strong, 
so that lepton asymmetries cannot increase A^cft signif- 
icantly above 3.4. Under those circumstances, if the cos- 
mological data (other than BBN) continues to push for 
large values of A^cff , new pieces of physics such as sterile 
neutrinos will be necessary to explain that excess. 
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